Nonlocal stabilization of nonlinear beams in a self-focusing atomic vapor 
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We show that ballistic transport of optically excited atoms in an atomic vapor provides a nonlocal 
nonlinearity which stabilizes the propagation of vortex beams and higher order modes in the presence 
of a self-focusing nonlinearity. Numerical experiments demonstrate stable propagation of lowest and 
higher order vortices over a hundred diffraction lengths, before dissipation leads to decay of these 
structures. 

PACS numbers: 42.65.Tg,42.65.Sf,32.80.-t 



The propagation and dynamics of localized nonlinear 
waves is a subject of great interest in a range of physi- 
cal settings stretching from nonlinear optics to plasmas 
and ultracold atomic gases [H, 0| . The structure and sta- 
bility of nonlinear optical modes is determined by the 
interplay of the radiation field with the functional form 
of the material nonlinearity In the case of opti- 

cal beams the nonlinear response can be described in 
terms of the induced change in the refractive index n 
which is often approximated as a local function of the 
wave intensity, i.e. n(r) = n(/(r)). However, in many 
real physical systems the nonlinear response is spatially 
nonlocal which means that the refractive index depends 
on the beam intensity in the neighborhood of each spa- 
tial point. This can be phenomenologically expressed as 
n(r) = J dr' K(r, r')J(r'), where the response kernel K 
depends on the particular model of nonlocality Q ■ 

It has been shown that nonlocality drastically affects 
the stationary structure and dynamics of spatial solitons, 
leading to such effects as collapse arrest of high intensity 
beams and stabilization of otherwise unstable complex 
solitonic structures @, 0,0, II- Nonlocality is often the 
consequence of transport processes which include atom 
or heat diffusion in atomic vapors plasma and 
thermal media or charge drift in photorefractive 

crystals [Hj]. In addition long range interactions are re- 
sponsible for a nonlocal response in liquid crystals [lH or 
dipolar Bose Einstein condensates [14j . 

Hot atomic vapors are an important and widely used 
nonlinear medium. The nonlocal nonlinear response of 
atomic vapors has previously only been associated with 
state dependent transport of ground state atoms which 
possess a multilevel structure [9| . In this letter we intro- 
duce a new and significant mechanism of nonlocality in 
atomic vapors which is provided by the ballistic trans- 
port of excited atoms and is important even for the sim- 
plest case of an idealized two-level atom. We show using 
parameters representative of beam propagation in Ru- 
bidium vapor that ballistic transport plays a dramatic 
role leading to stabilization of otherwise unstable vortex 
modes in the presence of a self-focusing nonlinearity. 



Prior to introducing a model for the nonlocal character 
of the refractive index we first recall the main features of 
beam propagation in a hot atomic vapor. We consider 
a scalar traveling wave E = e *(kz—wt) _|_ c c p or 

all parameters of interest the refractive index is n ~ 1 so 
the wave intensity is I ~ ^|£| 2 . In the slowly varying 
envelope approximation the paraxial wave equation is 
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where k — oj/c. The susceptibilities x' an d x" de- 
pend on atomic parameters. We assume a two-level 
atomic model for which the scattering cross section is 
a = (3A 2 /2tt) (1 + 4A 2 ) / 7 2 + I/I S )~ 1 , and the index 
of refraction is n — 1 — n a (a/k) (A0/7), where A Q is 
the transition wavelength, n a is the atomic density, 7 
is the full width at half maximum (FWHM) natural 
linewidth, Ao = uj — Lu a is the detuning between the 
optical frequency to and the atomic transition frequency 
Lu a = 27rc/A a , and I s is the saturation intensity. For a 
probe beam propagating along I in a hot vapor we make 
the replacement A — > A = A — kv z . Averaging over 
a Maxwell-Boltzmann velocity distribution at tempera- 
ture T gives an expression for the complex susceptibility 
which can be separated into a constant part and a part 
which depends on intensity. The results can be written 
in the form [1 51 ] x'q = xoInr[iv(a + ib)], 

XU 1 ) = Xo {Rc[Z(a + ibi)] - Rc[Z(a + ib)]} , (2a) 

^ = J^l-Mz(, + a)l|, (2b) 

[ V 7 ! + 1 1 Is J 

where xo = n a 67rbc 3 /io 3 , a = 2\/ln 2Ao/wd, b ■ 
\/ln 2j/ujd, bi — byl + I/I s , andwc = ky / 8\n2kBT/m 
is the FWHM of the Doppler profile for an atom of mass 
to. The plasma dispersion function is given by Z(z) = 
iy^e- z2 Erfc(~iz) where Erfc(z) = 1 - (2/0?) £ dt e"* 2 
is the complementary error function. 

For broad optical beams numerical solutions of Eqs. 
([I]) and ((2]) give an accurate description of propagation 
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FIG. 1: (color online) Characteristic length scales l\?~ 9 \ 
(,i 3 ~ e \ and Id in a Rb vapor cell as a function of tempera- 
ture. The inset shows the logarithm of the ballistic response 
function ~ (td/r) / °° d£ e~ r ^ ld ^e~^ (solid blue line) and 
the response function Ko(r/ld) for a 2D diffusive nonlocal 
equation [j| (dashed red line) as a function of r/ld- The re- 
sponse functions have been scaled to be equal at r — Id- 



effects in an atomic vapor. The physical effect leading to 
the nonlinear optical response is the transfer of popula- 
tion from the ground to the excited state and the creation 
of coherence between these states. Although motional ef- 
fects are accounted for as regards the Doppler smearing 
of the transition frequency, atomic motion also results 
in transport of excited atoms which leads to a nonlocal 
response. The degree of nonlocality depends on the char- 
acteristic length scales associated with the transport of 
excited state atoms. The first length scale is the mean 
free path atoms travel before a Rb-Rb collision occurs. 
This is given by £ c — (V^nacr)" 1 , where a is the col- 
lisional cross section. For collisions of ground state Rb 



atoms we use 16 ov_„ = 2.5 x 10 
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The cross 



section for collisions between excited and ground state 
collisions is much larger since these collisions occur via 
a long range dipole-dipole interaction [13]. The energy 
averaged cross section is <J g ~ e ~ 1,8> ^ — K 1 / 2 m 2 . The 
second length scale is Id = vt, which is the distance 
traveled by an atom moving at the most probable speed 
v = ^j2ksT /m in the 1 /e lifetime r of the excited state. 
For the 87 Rb 5P 3 / 2 level r ~ 26 ns. Figure [1] shows that 
for T < 155°C the ballistic transport length for excited 
atoms is l d ~ 7.5 fim and l d < l{ 9 ~ e] ,l [ ^ g) . Thus for 
these temperatures transport of excited atoms is ballis- 
tic with a length scale of Id- We note that the density 
at T = 155°C is n Q = 10 20 m -3 which is several orders 
of magnitude smaller than densities for which nonlocal 
effects due to the Lorentz local field are important [18j|. 

We wish to find an expression for the nonlocal material 
response that depends on the parameters r and v. We 
write the total atomic density n a = n g + n e as the sum of 



ground and excited state partial densities and introduce 
rate equations of the form 
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+ jn e + C g [n g ] 
jn e + C e [n e ], 
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where C g , C e are, as yet unknown, linear operators for 
ground and excited state atoms and a(T) = k[x'o+x"i{I)] 
is the absorption coefficient. If we assume that the total 
density is unchanged by the presence of the laser field 
(this is a reasonable assumption in hot vapors, but not 
for cold atoms) we must have £ 9 [n g ] = —C e [n e ]. If ex- 
cited state transport were a diffusive process we would 
have £ e [^e] = -D e V 2 n e and on dimensional grounds 

D e ~{vTf/T~iy T . 

The situation in the ballistic regime is different. The 
collisionless Boltzmann equation for the density of ex- 
cited atoms is ^f- = ^f- + v • V r n e . Working within 
the paraxial approximation we are interested in the two- 
dimensional problem where r = xx + yy and v = 
v x x + v y y. The Green function is found by solving 



+ v • V r n e = 5{t)n (r, v) 



(4) 



The time dependent solution of ((4]) is n e (r, v, t) — 
H(t)no(r — vt, v) where dH/dt — S(t). Now consider 
the solution for no(r,v) = <$(ro)/(v), with /(v) = 
{m/2irkBT)e~ v l v the two-dimensional thermal veloc- 
ity distribution. The solution is 



G(r,t;r ,t ) = 6{r - r - v(t - i ))- 
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The Green function for the spatial response is then 

,-V /v 
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This calculation neglects the fact that the excitation de- 
cays with rate 7. We can account for this by including a 
factor of e rt in the Green function to arrive at 



G r (r,t;r ,io;7) 
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The Green function is parameterized by the decay rate 
7. Since the rate of local excitation is al/hu the spatial 
distribution of n e is given by 



n e (r,t) = dt / dr G r (T,t;r ,to;j) 

J —OO J 

x^M J(r0)to) . 



(8) 
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The response to a temporally constant field which is a 
delta function in space, I = Iq8(tq), is, with Eq. |7J), 



n e (r) = 



HUJ TTVl 



d£ e 



-ir/{v(,) e -e 



(9) 



The Green function for the steady state spatial response 
which has units of s/m 2 is thus 



G(r,r ; 7 ) = / d£ e -W«) e -« 



(10) 



where r = |r — ro|. The result is plotted in the inset of 
Fig. [1] as a function of the scaled coordinate r/ld- We 
see, not unexpectedly, that the ballistic response falls off 
much more rapidly than the diffusive response. Note that 
J drG(r,0,7) = 1/7, since the time integrated response 
exponentially weights the input over a time window r = 
I/7. The spatial Fourier transform of the Green function 

is given by T[G] = ^ ^'^f Erfc[l/(fcl d )] which is well 
behaved with lim^^o F[G\ — 1-1 1- 

To complete the theoretical formulation of the wave 
propagation problem we need to calculate the nonlocal 
structure of the susceptibility Xni- The imaginary part 
of the nonlinear susceptibility is proportional to the dif- 
ferential density of excited and ground state atoms which 
decays with rate 7. When the intensity is uniform in 
space the susceptibility satisfies the relaxation equation 
dx'ni/dt = — 7(Xni — x"i) wnere the overbar denotes the 
steady state value of the susceptibility. When I is spa- 
tially varying we can use the Green function to write the 
stationary response as 



Xni( r )=7Xo / dr G(r,r ;7) 



\ v/l + 7(r )// s 



(11) 



The real part of the susceptibility is proportional to 
the coherence between ground and excited states which 
decays with rate 7/2. The Green function to be used for 
X^i is thus G(r,r ,7/2) and we can write 

Xni(r) = 1 f L J * G(r,r ; 7 /2) 

x {Re[Z(a + ibi[I{r )])] - Rc[Z(a + ib)}} . (12) 

Equations (flTl [T3J) together with the Green function Eq. 
dTUJ) and the wave equation |T]) are the main theoreti- 
cal result of this paper. They constitute a full descrip- 
tion of time-independent wave propagation in a two-level 
atomic vapor including Doppler broadening and trans- 
port induced nonlocality. 

The question of whether or not the ballistic transport 
is sufficient to stabilize nonlinear modes can be investi- 
gated by beam propagation calculations. We use param- 
eters corresponding to off-resonant propagation in a high 




FIG. 2: (color online) (a) Saturation function f sa t = He[Z(a+ 
ibi)] - Re[Z(a + ib)] for a = 4 and b = 0.0083. (b) Nonlocal 
single charged vortex mode with power P ~ 0.4 W (red line 
and red axis). The solid blue line shows the nonlocal nonlinear 
index Xni computed from Eq. (|12|l . the dashed blue line the 
local one computed from Eq. (|2a|l (blue axis). 



temperature Rb cell (A 
A /2tt = 1.46 GHz, I s 



10 



20 



780 nm, 7/2tt 
= 16.7 W/m 2 , 



= 6.07 MHz, 
T ~ 155°C, 



, id = 7.5 /im), which result in the dimen- 
sionless parameters a = 4, b — 0.0083, and xo = 0.03. We 
used this set throughout all the simulations presented in 
this paper. In the conservative system (xo = Xni = 0) 
all modes we tried (ground state, single charged vortex, 
dipole, double charged vortex) turned out to be stable if 
the power P = J I(r)dr, is high enough. At least for the 
latter two modes this is quite remarkable, since they are 
known to be unstable (or only stable in a small power 
window) for other nonlocal models d, [l9|, [2(|. We at- 
tribute this enhanced stabilization to the combination of 
nonlocality and nonlinear saturation. In fact, we inserted 
an artificial nonlinear saturation in the nonlocal thermal 
model used in Ref. [l9| and found that the double charged 
vortex becomes stable as well. However, it is worth point- 
ing out that nonlinear saturation without nonlocality does 
not stabilize higher order nonlinear modes [2l| . The local 
Eqs. (11I2|) feature a stable ground state only. 

Figure [5] illustrates both saturation and nonlocality for 
the nonlocal single charged vortex mode. If we consider 
only the saturation effect shown in Fig. [Ha) the result- 
ing nonlinear index is the dashed blue line in Fig. [2{b) • 
Together with the nonlocal kernel G [red line in Fig. [1] 
inset] we get the solid blue line, showing some filling in 
of the central dip in the index profile, and the formation 
of a broader "nonlocal waveguide" . 

The limiting mechanism with respect to long distance 
propagation of higher order nonlocal nonlinear modes is 
not destabilization but dissipation. The action of both x'o 
and y", is not negligible over one diffraction length Zd — 
2ki d [22J. As an illustrative example, the propagation 
of the nonlocal single charged vortex mode, is shown in 
Fig. EJa). As input power we use about 0.4 W. Note 
the clearly visible influence of the nonlinear term Xni in 
the blue power curve. The nonlocal vortex survives a 
propagation distance of more than 150zd [see Fig. EJb)]. 
For comparison, the propagation of the local vortex with 
the same input power is shown in dashed lines in Fig. 
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FIG. 3: (color online) (a) Nonlocal (solid lines) and local 
(dashed lines) dissipative propagation of the single charged 
vortex mode with input power 0.4 W. The blue lines and 
blue axis show the beam power, the red lines and red axis the 
maximal intensity versus propagation distance, (b) Intensity 
and phase distribution of the nonlocal single charged vortex 
at input z — and at z = 160zd just before it decays, (c) 
Maximal intensity (red) and FWHM (green) of the nonlocal 
single charged vortex as a function of beam power. The solid 
lines are computed upon propagation, the diamonds from sta- 
tionary numerical solutions of the conservative problem, (d) 
Intensity and phase distribution of the local single charged 
vortex at z = and at z — 15zd when it decays, (e) Intensity 
and phase distribution of the nonlocal dipole mode at z = 
and at z = 160zd just before it decays, (f) Same for the 
nonlocal double charged vortex at z — and at z — WOzd- 



OH a). This vortex disintegrates after less than 15 Zd [see 
Fig. 02d)]. Hence, we clearly see that the stabilization is 
due to nonlocality. With the same input power of about 
0.4 W, we also observed a robust nonlocal dipole [see Fig. 
EHe)] and double charged vortex [see Fig. [3£f)]. 

The key feature enabling robust nonlocal dissipative 
propagation over a hundred diffraction lengths is the 
above mentioned stability for high powers. Starting in 
the stable power regime, dissipation makes the nonlin- 
ear mode "glide down" the family branch until it reaches 
powers in the unstable regime. Figure 03c) confirms this 
property by comparing maximal intensity and FWHM 
obtained upon propagation with values found from ex- 
act numerical solution of the conservative problem using 
the method described in 0] . The solid lines are obtained 
upon propagation, which explains the small oscillations 
in the curves. 

In conclusion, we have shown that ballistic transport 
of optically excited atoms in a thermal vapor provides 



a generic nonlocal nonlinearity which can stabilize the 
propagation of vortices and other higher order modes in 
a self-focusing medium. For sufficiently high power we 
found a stable dipole mode and single and double charged 
vortices. In realistic models dissipation is not negligible. 
Nevertheless, numerical experiments demonstrate robust 
propagation over a hundred or more diffraction lengths. 
This is possible due to adiabatic conversion into solitons 
with lower power, but of the same family. 

Numerical simulations were performed on the SGI Al- 
tix 3700 Bx2 cluster of the Australian Partnership for Ad- 
vanced Computing and on the IBM p690 cluster (JUMP) 
of the Forschungs-Zentrum in Julich, Germany. 
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